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We report the existence of metallic two dimensional electron gases (2DEGs) at the (001) and (101) 
surfaces of bulk-insulating Ti 02 anatase due to local chemical doping by oxygen vacancies in the 
near-surface region. Using angle-resolved photoemission spectroscopy, we find that the electronic 
structure at both surfaces is composed of two occupied subbands of dxy orbital character. While 
the Fermi surface observed at the (001) termination is isotropic, the 2DEG at the (101) termination 
is anisotropic and shows a charge carrier density three times larger than at the (001) surface. 

Moreover, we demonstrate that intense UV synchrotron radiation can alter the electronic structure 
and stoichiometry of the surface up to the complete disappearance of the 2DEG. These results open a 
route for the nano-engineering of confined electronic states, the control of their metallic or insulating 
nature, and the tailoring of their microscopic symmetry, using UV illumination at different surfaces 
of anatase. 


In its pure stoichiometric form, the transition metal 
oxide (TMO) Ti02 is a transparent insulator that crys¬ 
tallizes in mainly two different phases: rutile and anatase. 
Both phases have been studied extensively over the last 
decades, due to their photocatalytic properties discussed 
in several reviews W- Recently, a strong interest in 
the anatase phase of Ti 02 also surged, owing to its po¬ 
tential for applications in other research fields. For in¬ 
stance, networks of anatase nanoparticles are found in 
dye-sensitized solar cells m, anatase thin films can be 
used as transparent conducting oxides [7], and devices 
based on anatase can be envisioned in spintronics 0111. 
To harness such a wide range of functionalities and guide 
potential applications using anatase, it is thus critical 
to understand its microscopic electronic structure, which 
will be ultimately responsible for the remarkable proper¬ 
ties of this material. Moreover, as most applications in 
microelectronics or heterogeneous catalysis involve essen¬ 
tially the electronic states at the material’s surface, it is 
crucial to directly measure and characterize such states. 

More generally, the study of two-dimensional electron 
gases (2DEGs) in TMO surfaces/interfaces has become 
a very active field of research. The archetypal example, 
SrTiOs-based heterostructures, display many fundamen¬ 
tally interesting properties [lOHig, such as field-effect in¬ 
duced insulator-to-superconductor transitions US], mag- 
netism [14] and the coexistence of magnetism and su¬ 
perconductivity m- More recently, the discoveries of 
2DEGs at the bare (001), (110) and (111) surfaces of 
SrTiOs [T 6 II 2 T] and KTaOs [221(24] . triggered new av¬ 
enues of research. Additionally, very recent spin-resolved 
measurements of the 2DEG on SrTiOs (001) revealed a 
giant spin splitting of its electronic structure [25|, epit¬ 
omizing the richness of physical properties that can be 
found in these 2DEGs. 

It is well established that the TiO atomic planes, and 


their ability to accommodate chemical doping by oxygen 
vacancies at the surface region, play a key role in the 
formation of the 2DEG at the surface of SrTiOs (001). 
Thus, as a step forward to understand the formation of 
2DEGs in TMOs, it is natural to focus on pure Ti02 crys¬ 
tals such as rutile or anatase. In fact, it is known that, 
for both rutile and anatase crystal surfaces, UV or elec¬ 
tron irradiation creates oxygen vacancies [261(28] . But 
while 2DEGs can be readily obtained at various anatase 
surfaces, as we will show further, the failure to observe 
2DEGs on the cleaved or in-situ prepared surfaces of ru¬ 
tile Ti02 (110) (Ref. [29] and Supplemental Material), 
indicates that other factors are equally important. In 
fact, excess electrons in Ti 02 , due to oxygen vacancies 
at the surface, form polarons which behave quite differ¬ 
ently for the rutile and anatase phases m- This demon¬ 
strates that structural factors, such as surface lattice dis- 
torsions, can determine the fate of those electrons. Thus, 
to elucidate the origin and realization of 2DEGs in differ¬ 
ent oxides, it is essential to probe the surface electronic 
structure arising in different lattice configurations. 

Here we report the observation of 2DEGs at the (001) 
and (101) surfaces of Ti02 anatase. Using angle resolved 
photoemission spectroscopy (ARPES), we find that, for 
both surfaces, the electronic structure consists of two 
light bands of dxy character. While the Fermi surface 
at the (001) surface is isotropic, the one at the (101) sur¬ 
face is highly anisotropic, reflecting the symmetry of the 
surface layer m- Moreover, we show that the character 
of the excess electrons (localized or delocalized) depends 
on the concentration of oxygen vacancies, and that the 
2DEG disappears, leaving no mobile conduction electrons 
near the Fermi level, if the surface is reduced too much 
by intense UV irradiation. 

The ARPES measurements were conducted at the Syn¬ 
chrotron Radiation Genter (SRG, University of Wiscon- 
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FIG. 1. (Color online) Angle integrated 
spectra of (a) the Ti-3p peak, and (b) the 
in-gap state of the anatase Ti 02 (001) and 
(101) surfaces (upper and lower curves, re¬ 
spectively) measured at hu = 100 eV. The 
black curve was measured shortly after the 
first light exposure, and the blue, red and 
green curves at sequentially increasing times, 
(c) Binding energy of the Ti^+ peak and 
concentration of oxygen vacancies (ratio be¬ 
tween the Ti^+ and total Ti-3p peak area) 
for different durations of irradiation expo¬ 
sure of the (101) surface (upper panel), and 
binding energy and intensity of the in-gap 
state (lower panel). Analogous results, not 
shown, are obtained for the (001) surface 


sin, Madison) and the CASSIOPEE beamline of Syn¬ 
chrotron Soleil (Erance) on surfaces cleaved at low tem¬ 
peratures under ultra-high vacuum - see Supplemental 
Material for details about the sample preparation and 
measurements. All through this paper, directions and 
planes are defined in the simple tetragonal conventional 
cell of anatase. We note [hkl] the crystallographic direc¬ 
tions in real space, (hkl) the corresponding directions in 
reciprocal space, and (hkl) the planes orthogonal to those 
directions. The indices k, and I of T^ki correspond to 
the reciprocal lattice vectors of the body-centered unit 
cell. 

As shown in Eigs. [^a) and Bb), the presence of oxy¬ 
gen vacancies in the surface region of anatase can be iden¬ 
tified in the photoemission spectra by the formation of a 
shoulder in the Ti-3p peak related to a lower valence of 
the Ti cation, and the appearance of an in-gap state m 
around 1 eV below the Eermi level corresponding to 

electrons trapped near oxygen vacancies [30], The black 
curves were measured shortly after the first exposure of 
the sample to the synchrotron radiation, while the blue, 
red and green curves were recorded at later subsequent 
times specified respectively by the abscissas of the blue, 
red and green open squares in Eig. [^c). 

As seen in Eigj^a), the peak position of the Ti-3p peak 
shifts to higher binding energies upon UV irradiation, 
demonstrating the band bending (bb) at the surface. We 
fit the Ti-3p peak by one or two Voigt functions plus 
a Shirley background to determine the contribution of 
Ti^+ and Ti^+ states to the line shape. This yields a 
shift of Ti^+ peak of = 110 meV for the (001) 

surface, and = 150 meV for the (101) surface, 

which corresponds to the minimal band bending at the 
surfaces averaged over the photoemission probing depth 
at the used photon energy. The actual band bending 
might be larger, as oxygen vacancies might be already 


induced at the surface by the cleaving process or by the 
short beam exposure before the first measurement. 

The red and green curves in Eigs. [^a) and[^b), show 
that many more oxygen vacancies can be induced at the 
Ti 02 (101) surface than at the (001) surface, in agree¬ 
ment with previous studies [32] . The band bending (peak 
position of the Ti^+ peak) and concentration of oxygen 
vacancies (ratio between the Ti^+ and total Ti-3p peak 
area) for the blue, green, and red curves as well as for 
intermediate measurements are plotted in figure [^c) for 
the specific case of the anatase (101) surface. This shows 
that the band bending saturates at rather low concentra¬ 
tion of vacancies at the surface. In fact, recent theoretical 
studies on the electronic structure of SrTiOs [33] [34] sug¬ 
gest that, for a high concentration of oxygen vacancies, 
excess electrons do not fill up the conduction band any 
longer, but form only localized states, and thus do not 
contribute to the rigid band shift. Additionally, as we 
will see later, the 2DEG disappears if a high concentra¬ 
tion of oxygen vacancies is created. 

Similarly, as shown in figures Bb) and (c) (lower 
panel), the intensity and binding energy of the in-gap 
state evolve significantly upon UV irradiation, both for 
the (001) and (101) surfaces. Possible explanations for 
the different binding energies are that the oxygen vacan¬ 
cies are located at different lattice sites and/or forming 
clusters [33], or the existence of different types of de¬ 
fects [351 [36] . 

Thereafter, we focus on the electronic structure of the 
2DEG at the Eermi level, obtained once a stable band¬ 
bending has been attained after weak exposure to UV 
irradiation (typically ^ 20 min) - see Eig. Be). Fig- 
ure|^a) shows the bulk Brillouin zone of anatase, its high 
symmetry directions, and the planes of measurements at 
its cleaved (001) surface. The Eermi surface at the (001) 
plane, shown in figure j^b), consists of two circles of dxy 
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FIG. 2. (Color online) (a) Bulk Brillouin zone of anatase, showing the relevant high symmetry directions (green) and planes (red 
and blue) of the Fermi surfaces and band dispersion at the (001) surface, (b) Fermi surface map (second derivative of ARPES intensity, 
negative values) measured at hu = 47 eV on a cleaved, bulk insulating Ti02 anatase (001) surface. The red lines indicate the edges of 
the bulk Brillouin zones in the (001) plane around roi2- The same relative color scale indicating minimum (min) to maximum (max) 
ARPES intensities is used in all color figures in the rest of this work, (c) Fermi surface map (second derivative of ARPES intensity, 
negative values) in the fc^ooi) “ ^( 010)5 (100) plane, acquired by varying the photon energy in 1 eV steps between hui = 38 eV and 

hh'2 = 85 eV. To calculate the momentum perpendicular to the surface we set the inner potential to Vq = 13 eV EH 133 . The blue lines 
are the bulk Brillouin zones containing rio2 and Pios- (d) Energy-momentum map measured at hu = 47 eV along the (010) direction on 
the (001) surface. Dashed red lines are parabolic fits to the dispersing light bands. Red arrows show a kink attributed to electron-phonon 
coupling |29| . (e-h) Same as (a-d) on the anatase (101) surface. Data in panels (f) and (h) were measured at hiy = 47 eV. The green 
dashed arrow in (f) indicates the direction of measurements along ky in (g) and (h), which is slightly off (010). The out-of-plane Fermi 
surface map (second derivative of ARPES intensity, negative values) in (g) was measured by varying the photon energy in 1 eV steps 
between hui = 33 eV and hh '2 = 51 eV and between hus = 76 eV and hv '4 = 103 eV. The Supplemental Material includes figures for the 
raw spectra, Fermi surface maps, and the periodicity of the electronic structure in the (001) and (101) planes. 


orbital character, confirmed by light polarization depen¬ 
dent measurements presented in the Supplemental Mate¬ 
rial, similar to the case of SrTiOs p^UTS] . The quasi-2D 
character of these two electronic states is strictly demon¬ 
strated by the Fermi surface map in the (010) — (001) 
plane, perpendicular to the surface plane, shown in fig¬ 
ure [^c). This Fermi surface shows that the bands are es¬ 
sentially non dispersing along /c^ooi) over a bulk Brillouin 
zone, thereby confirming the confined character of the 
electrons close to the surface. Figure [^d) presents the 
energy-momentum map close to the bulk F 102 point along 
the cut parallel to the (010) direction shown in panel (a). 
Two dispersive light bands are visible. Their bottoms 
are at —60 meV for the upper band, and —172 meV 
for the lower band. The Fermi momenta of 0.09 
and 0.16 and the experimental band dispersions are 
given by the peak positions of the momentum distribu¬ 
tion curves (MDGs) at and below Ep. A parabolic fit to 
the band dispersion yields an effective mass of approxi¬ 
mately 0.5me, which agrees well with theoretical [38l [39] 
and experimental [40] results on bulk Nb-doped anatase. 


Assuming two spin orientations per band, this gives an 
electron concentration of ^ 5.4 x 10^^ cm“^, or 

about 0.08 electrons per a^, where a is the short lattice 
constant of the tetragonal lattice and the cross section 
of the unit cell in the (001) plane. Note that a kink in the 
dispersion at a binding energy of about 70 meV is visible 
(red arrows). In the case of bulk states of anatase, such 
a kink has been attributed to electron-phonon coupling 
and studied in detail by Moser et al. [29] . 

Next, we present the ARPES measurements on the 
(101) surface of anatase. The bulk Brillouin zone is 
shown again for clarity in Fig.j^e), together with relevant 
directions for the (101) surface and the planes of the mea¬ 
sured Fermi surfaces. Figure [^f) shows the Fermi surface 
on the (101) plane. It consists of two ellipses of identi¬ 
cal shape, corresponding again to d^y orbitals, similar to 
the dxy derived ellipses observed for the 2DEG at the 
SrTiOs (110) surface [T9l[20]. The Eermi surface map in 
the ky — kz plane, perpendicular to the cleaved surface, 
is shown in figure |^g) . The bands are essentially not 





















4 


dispersing along kz^ indicating their 2D character. Note 
that, in contrast to the (001) surface, the intensity of the 
2DEG states in the kz — ky plane at the (101) surface 
drops quickly far from bulk F points, due to final-state 
effects in the photoemission cross-section (see the discus¬ 
sion in the Supplemental Material of [20]). Figure [^h) 
shows the energy-momentum map close to the bulk r 2 io 
point along the ky direction. As in the (001) surface, 
two light bands, upper and lower, are also observed in 
this case. Their bottoms of are located at —60 meV and 
— 130 meV, respectively. The Fermi momenta taken from 
the Fermi surface in Fig. [^f) are 0.41 & 0.33 

in the (101) direction and 0.15 A“^ & 0.10 A“^ in the 
(010) . The parabolic dispersion based on the experi¬ 
mental Fermi momenta and band bottom energies yields 
an effective mass of ^<oio> = O.GrUe along the < 010 > 
direction, very close to the one measured at the (001) 
surface along the same direction. Assuming again 2 
spins per band, this gives an electron concentration of 
~ 1.5 X 10^^ cm“^, three times larger than the 
one at the (001) surface, which is probably related to the 
higher concentration of oxygen vacancies at the (101) sur¬ 
face [32] . 

The presence of two non-degenerate bands of identical 
orbital character directly implies that the probed elec¬ 
tronic structure is not simply the one expected for bulk 
anatase. For SrTiOs, the existence of different subbands 
has been related to the confinement of the electronic 
states in a quantum well at the surface [laiTT]. More 
recently, spin-resolved ARPES experiments related the 
two dxy-hands in SrTiOs to different spin polarizations 
of the ground-state dxy subband [25]. Whether or not a 
magnetic order is also present at the oxygen-deficient sur¬ 
face of Ti02 anatase cannot be deduced from the ARPES 
data presented in this paper. Note that magnetic order 
in thin films of Ti 02 due to oxygen vacancies was previ¬ 
ously observed and discussed [4T]-[46], and is thus a po¬ 
tential reason for the observed splitting of the da^^-bands. 
Further experimental work is necessary to clarify the pos¬ 
sible magnetism at the oxygen-deficient surface of Ti 02 
anatase. 

In contrast to the 2DEGs at the surface of SrTiOs and 
KTaOs, OUT ARPES data in Ti02 anatase do not show 
any heavy {dzx or dyz) subbands -which we confirmed 
by complementary measurements as a function of pho¬ 
ton polarization and energy, shown in the Supplemental 
Material. The degeneracy of the three t 2 g orbitals is lifted 
in bulk anatase m due to its tetragonal crystal struc¬ 
ture and the inequivalent Ti-0 bonding in the x — y plane 
compared to the 2 ; direction. Gonsequently, the electron 
gases confined at the surfaces of Ti 02 anatase are only 
composed of electrons of dxy-OTh\tdl character. As shown 
in the Supplemental Material, the spatial extent of the 
2DEGs in anatase is similar to the ones at the surfaces 
of SrTiOs and KTaOs, namely about 2 nm. 



FIG. 3. (Color online) Energy-momentum maps measured at 
hi/ = 47 eV along the (010) direction at the anatase (001) surface 
for different UV doses, (a) After a few minutes of UV irradiation, 
(b) After several hours of exposure to UV light, (c) After moving 
the UV spot to a region in the sample not irradiated before, and 
illuminating again for a few minutes. The white curves correspond 
to the MDCs integrated over Ep ± 10 meV, and are plotted using 
the same absolute vertical scale in the three panels. 


We finally show that exposure to a large UV irradi¬ 
ation dose leads to the local destruction of the 2DEG. 
Figure [^a) displays the energy-momentum map around 
F 012 measured in a freshly-cleaved anatase (001) surface 
after only a few minutes of UV irradiation. The 2DEG 
is clearly visible. However, as shown in Fig. [^b), the 
spectral weight of the dispersing features essentially dis¬ 
appears after 1-2 hours of exposure to UV light, leaving 
almost no intensity at the Fermi level. As seen in fig¬ 
ure [^c) , the 2DEG is fully recovered after moving the 
UV spot to a neighboring region in the sample, not ir¬ 
radiated before, and illuminating for a few minutes. As 
shown in the Supplemental Material, a high concentra¬ 
tion of oxygen vacancies induced by the UV irradiation 
also destroys the crystallinity at the surface, which might 
in turn lead to the localization of the itinerant electrons 
forming the 2DEG. As the size of the beam spot is only 
~ 50 X 50 /im^, the strong reduction of the surface is 
local. 

In conclusion, we demonstrated the existence of 
2DEGs at the (001) and (101) surface of Ti02 anatase. 
The two dxy subbands composing the 2DEGs form con¬ 
centric circular Fermi surfaces at the (001) termination 
and ellipsoidal Fermi surfaces at the (101) termination. 
Such orientational tuning of the Fermi sea symmetries is 
wholly analogous to the one found at different surfaces 
of SrTiOs [20]. Furthermore, we found that light irradia¬ 
tion locally dopes the anatase surfaces, while the behav¬ 
ior of the excess electrons, i.e. delocalized or localized, 
depends on the concentration of oxygen vacancies at the 
surface, eventually resulting in the local destruction of 
the 2DEG after a high UV irradiation dose. This effect, 
never observed for 2DEGs at other oxide surfaces, shows 
that UV light can be used to tailor locally the surface or¬ 
der/disorder in anatase, and thus engineer nano-patches 
of metallic 2DEG alongside patches of disordered insu¬ 
lating material. 
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SUPPLEMENTAL MATERIAL 


Methods for sample preparation and measurements 

The non-doped, polished crystals of Ti02, of typical 
sizes 5 X 5 X 0.5 mm^, were supplied by SurfaceNet GmbH. 
The anatase crystals were natural grown and of orange 
color, with less than 5 ppm of Mn impurities. The orien¬ 
tation of the large surface of the crystals was either (101) 
(natural cleaving plane) or (001) for anatase, and (110) 
for rutile. Clean and crystalline surfaces were obtained 
by cleaving the crystals in-situ^ using the standard top- 
post procedure m, at temperatures between 7 K and 
25 K and pressure lower than 6 x 10“^^ Torr. To provide 
electrical grounding to the exposed surface, the crystals 
were glued with a conducting silver epoxy to the sample 
holder and, after gluing the top-post, wholly covered with 
graphite paint. The same protocol has been successfully 
used in previous ARPES studies of other oxides, such 
as SrTiOs and KTaOs, that are transparent and highly 
insulating in the bulk ng ng. The sharp photoemis¬ 
sion lines and dispersing bands in the ARPES measure¬ 
ments demonstrate a crystalline surface. For the different 
cleaved surfaces, the periodicity of the measured elec¬ 
tronic structure, discussed later (figure [^, corresponds 
to that of an unreconstructed surface, and confirms the 
surface orientation. 

The ARPES measurements were conducted at the Syn¬ 
chrotron Radiation Center (SRC, University of Wiscon¬ 
sin, Madison, USA) and the CASSIOPEE beamline of 
Synchrotron SOLEIL (Saint-Aubin, France) We used lin¬ 
early polarized photons in the energy range 38 — 103 eV 
and hemispherical electron analyzers with vertical slits. 


The total angle and energy resolutions were 0.25° and 
25 meV at SRC, and 0.25° and 15 meV at Soleil. The 
mean diameter of the incident photon beam was smaller 
than 50 fim at SOLEIL, and about 200 fim at SRC. The 
photon flux at SOLEIL was ^ 2 x 10^^ photons/s/0.Ibw, 
and about 3 times smaller at SRC. Thus, the density 
of photons irradiating the sample was approximately 50 
times larger at SOLEIL than at SRC. The results were 
reproduced for at least three different samples for each 
surface orientation. 


Rutile vs. anatase 

The different behavior of excess electrons induced by 
oxygen vacancies in rutile and anatase was already dis¬ 
cussed in the main text. In contrast to anatase, the 
excess electrons in rutile do not fill up the conduction 
band (forming the 2DEG), but are localized at lattice 
Ti sites m- Accordingly, we did not detect any dis¬ 
persing features at the Fermi level at the surface of ru¬ 
tile (110), neither for cleaved nor for in situ prepared 
surfaces under various measurement conditions (photon 
energy, emission angle, light polarization) at T = 7 K. 
To illustrate the difference between the two systems, fig¬ 
ure |^a,b) shows the valence band, the in-gap state cor¬ 
responding to the localized electrons, and in the case 
of anatase the 2DEG measured at normal emission at 
hn = 100 eV. The angle-integrated data in figure [^c) 
demonstrates that the concentration of oxygen vacancies 
is similar in the two cases, as the in-gap states have simi¬ 
lar intensities. The failure to observe dispersing features 
at the Fermi level of rutile shows that oxygen vacancies 
at the surface are not sufficient to create a 2DEG, and 
that structural factors play a crucial role regarding the 
localization/delocalization of excess electrons. 

Surface vs. bulk-like electronic structure of anatase 

Moser et al [29] characterized the polaronic conduc¬ 
tion for different charge carrier densities (controlled by 
the oxygen partial pressure sample chamber) at the 
anatase TiO2(001) surface using ARPES. In the present 
work, we only show data after reaching the saturation 
value of the charge carrier density. Intermediate val¬ 
ues of the charge carrier density can be observed but 
are not stable long enough to measure with a sufficient 
signal-to-noise ratio without controlling the oxygen par¬ 
tial pressure in the measurement chamber. Contrary to 
our results, Moser et al observed only one dispersing 
subband of 3D character (z.e., dispersing along the di¬ 
rection perpendicular to the surface) at low charge car¬ 
rier densities. Note that for a low concentration of oxy¬ 
gen vacancies and charge carriers, the potential well at 
the surface is rather shallow and the electrons are thus 
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FIG. 4. (Color online) (a,b) Energy-momentum maps measured at hu = 100 eV along the (010) or the (001) direction at the anatase 
(101) and rutile (110) surface. The intensities in different energy intervals (valence band, in-gap, 2DEG) were normalized to make all the 
features visible, (c) Angle-integrated, raw spectra of the two energy-momentum maps in (a,b). The inset shows a zoom of the angle- 
integrated spectra (log scale) for angles close to normal emission. The spectral weight at the Fermi level in the case of the anatase (101) 
surface corresponds to the 2DEG. 


barely confined. But, even for charge carrier densities 
comparable to our results {tisd,M oser ~ 2 x 10^^ cm“^ 
to nsD ~ ~ 4 X 10^^ cm“^) Moser et al 

observed only one dispersing band. Such discrepancy 
in the number of bands and dimensionality occurred al¬ 
ready in the case of the electronic structure at the sur¬ 
face of SrTiOs: Chang et al [48], using the same ARPES 
setup as Moser et al, measured only one dispersing band 
of dxy character and attributed it to the bulk, whereas 
two dxy bands and 2D behavior were observed in other 
works [TBHTS] . Both studies observing only one dxy band 
were conducted at the same synchrotron endstation at 
the Advanced Light Source. Thus, the discrepancy in 
the number of bands might be related to the specifics 
of the ARPES endstation. Possible reasons are a lower 
photon flux or different photoemission selection rules due 
to different measurement geometries. Eor instance, for 
SrTiOs, we checked (not shown) that using horizontal 
slits instead of vertical slits for the electron analyzer 
yields only one of the two dxy subbands. 

Spatial extent of the 2DEG at the surface of anatase 

An estimate of the spatial extent of the 2DEGs in 
anatase using the triangular potential well model used 
previously [T6| can be obtained, subject to one caveat. 
If the subband splitting is due to the quantum confine¬ 
ment, the width of the potential well follows directly 
from this splitting and the effective mass of the sub¬ 
bands along the confinement direction m- Eor anatase, 
= 4.65me [39], yielding a confinement length of 
1.8 nm for the 2DEG at the anatase (001) surface. If, 
on the other hand, the subband splitting is due to mag¬ 
netism m, one can assume that the next quantum-well 


states of dxy character lie just above Ep, which gives an 
upper bound of 1.4 nm for the width of the potential well 
at the (001) surface. One sees that, in any case, the con¬ 
finement length of the 2DEGs at the surface of anatase is 
similar to the ones at the surfaces of SrTiOs and KTaOs, 
of about 2 nm. 


Influence of UV synchrotron irradiation 

We present now supplementary data regarding the in¬ 
fluence of the UV synchrotron irradiation on the sto¬ 
ichiometry and electronic structure at the surfaces of 
anatase. 

We conducted our measurements at two different syn¬ 
chrotrons: SRC {2^^ generation) and Soleil (3^ genera¬ 
tion). The normal measurement conditions (photon flux 
density at the surface of the sample, vacuum pressure, 
temperature) at SRC enable to reduce the surface, cre¬ 
ate the 2DEG, and study a stable 2DEG for days. In 
contrast, the much higher photon flux density at Soleil 
creates so many oxygen vacancies that the induced disor¬ 
der at the surface lead to the disappearance of the 2DEG 
within about one hour at T = 7 K at the anatase (101) 
surface. As this effect was never observed for the 2DEG 
at the surface of SrTiOs and KTaOs, we conclude that 
the anatase surfaces are much more sensible to the UV 
irradiation. 

The effects of the synchrotron UV light on the (001) 
surface are shown in figure An approximate upper 
limit for the effects is shown in Eig. [^a), where we 
present the angle-integrated spectra measured at hi' = 
100 eV on a fractured (001) anatase sample for bind¬ 
ing energies between —45 eV and 2 eV. The red curve 
corresponds to a weakly reduced surface, similar to the 
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FIG. 5. (Color online) (a) Angle-integrated spectra of a cleaved TiO 2 ( 001 ) anatase sample measured at hu = 100 eV showing the 
density of states for binding energies between —45 eV and 2 eV. The red curve corresponds to a weakly reduced surface, the blue one to 
a surface exposed to the intense zero-order light irradiation of the synchrotron, (b, c) Corresponding angle-resolved energy-momentum 
maps (second energy derivatives, negative values only) zoomed over the valence-band region. While the valence band is dispersive in 
the crystalline, weakly irradiated surface, it becomes non-dispersive and featureless in the amorphous, strongly irradiated surface, (d) 
Angle-integrated spectra measured at hu = 47 eV (log scale for spectral intensity), showing the tail of the in-gap peak and the intensity 
at the Fermi level. The red curve (label 1) corresponds to a freshly cleaved sample irradiated for a first sweep, the blue curve (label 2) 
to a surface exposed to irradiation for several hours, and the black curve (label 3) to the electronic structure after changing the sample 
position, (e, f, g) Energy-momentum maps measured at hu = 47 eV along the (010) direction. The momentum-integrated curves of these 
maps corresponds to the data shown in (b) as indicated by the labels 1, 2 and 3. 


data in figure 1 in the main text, and the blue curve 
to the electronic structure of the surface after exposing 
the sample to high intensity zero-order irradiation of the 
synchrotron. As can be seen from figure there are sev¬ 
eral changes in the electronic structure induced by the 
severe exposure to UV light. First, the shoulder of the 
Ti-3p peak, corresponding to Ti^+ states, transforms to 
become the the main Ti-3p peak at even lower binding 
energies, indicating that the majority of Ti atoms close 
to the surface have a low valency (< 3+) due to the cre¬ 
ation of many oxygen vacancies. Second, the shape of 
the valence band becomes featureless and its dispersion 
disappears after a high UV dose, demonstrating the large 
disorder induced at the surface. Third, a high concentra¬ 
tion of oxygen vacancies is obvious from the increased 
intensity of the in-gap peak and the decreased intensity 
of the peaks related to oxygen (0-2s and valence band). 


These measurements clearly demonstrate that the sur¬ 
face structure can be altered drastically, to the point of 
completely destroying the crystalline order, due to light 
irradiation. 

Figures l^b, c) show the angle-resolved measurements 
corresponding to the spectra in figure |^a), zooming over 
the valence-band region. The clear dispersion of the va¬ 
lence band in the weakly reduced surface. Fig. |^b), is 
a direct consequence of the good surface crystallinity 
(Bloch theorem). After intense UV irradiation, Fig.j^c), 
the surface crystallinity is destroyed, and the valence- 
band becomes non-dispersive and featureless. 

Even without the zero-order irradiation, the reduction 
of the surface can lead to the disappearance of the 2DEG. 
The angle-integrated spectra in figure |^b) were mea¬ 
sured at hu = 47 eV and show the tail of the in-gap 
peak and the intensity at the Eermi level for binding en- 
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ergies between —0.4 eV and 0.1 eV. After several hours 
of irradiation (blue curve), the intensity of the in gap 
state increases and the intensity at the Fermi level corre¬ 
sponding to the 2DEG decreases. The state of the weakly 
irradiated surface (red curve) can be retrieved by slightly 
changing the position of the UV light spot on the sam¬ 
ple surface (black curve). As the size of the beam spot is 
only ^ 50x50 , the strong reduction of the surface is 

local. This results suggests that UV light can be used to 
tailor locally the surface order/disorder in anatase, and 
thus write mesoscopic or nanoscopic patches of metallic 
2DEG intercalated with patches of disordered insulating 
material. 

The energy momentum maps corresponding to the 
angle-integrated spectra of figure ib) were already 
shown and discussed in figure 3 of the main text. We 
display them again, for completeness, in figures [^c,d,e). 
The 2DEG is clearly visible in figure [^c), whereas there 
is almost no intensity at the Eermi level in figure [^d) at 
the strongly reduced surface, and the 2DEG is recovered 
in figure HJ; e), after moving the UV spot to a neighboring 
region in the sample. 

Another interesting fact is that the anatase (101) sur¬ 
face degrades more rapidly than the (001) surface. The 
data corresponding to this degradation are shown in 
figure The angle-integrated spectra (red curves) in 
Eig. [^a) were measured at equal time intervals (of ap¬ 
proximately 2.5 minutes) at a photon energy hu = 47 eV, 
and show the increase of the intensity of the in-gap peak 
and decrease in the intensity at the Eermi level for bind¬ 
ing energies between —2.5 eV and 0.2 eV. Three of the 
curves highlighted in black, green and blue correspond to 
the start, a midpoint and the end of these measurements. 
The intensity of the 2DEG is evident in the momentum 
distribution curves at the Eermi level in figure |^b). At 
the beginning (black curve) two peaks, corresponding to 
the dispersing bands shown in figure 2(h) of the main 
text, are visible. Their intensity decays rapidly after sev¬ 
eral minutes, is barely visible (green curve) and disap¬ 
pears completely (blue curve). Note that the Eermi mo¬ 
menta in the MDGs does not increase -compare black 
and green curve in Eig. |^b), although more oxygen va¬ 
cancies are created as evidenced by the increase in in¬ 
tensity of the in-gap state in Eig. [^a). As the electron 
density in the conduction band is proportional to the 
Eermi momenta n oc kp, the excess electrons due to oxy¬ 
gen vacancies seem to populate the in-gap state (localized 
electrons) and not any longer the conduction band (delo¬ 
calized electrons). This observation is in agreement with 
the observed saturation of the band bending in the main 
text. 

A high concentration of oxygen vacancies induced by 
the synchrotron irradiation creates a high degree of disor¬ 
der at the surface. If the degree of disorder at the surface 
is sufficiently high, no dispersing bands exist as there is 
no longer the periodicity of the lattice, and the 2DEG 



FIG. 6. (Color online) (a) Angle-integrated spectra (log scale 
for intensities) measured at equal time intervals at hu = 47 eV for 
binding energies between —2.5 eV and 0.2 eV. One observes the in¬ 
gap peak and the intensity of the 2DEG at the Fermi level. Three 
of the curves are highlighted in black, green and blue corresponding 
to the start, a midpoint and the end of these measurements, (b) 
Corresponding momentum distribution curves at the Fermi level. 


disappears, as effectively seen in Fig. [^b). 


Raw data and periodicity of the electronic structure 

Figure shows the raw Fermi surface maps and en¬ 
ergy distribution curves (EDCs) of the data presented 
in Fig. 2 of the main text. Gomparing the raw maps 
of Figs, [^a, b, d, f) with the second-derivative maps of 
Figs. 2(b, c, f, g) in the main text, it is evident that 
the second derivative in the Fermi surface maps only en¬ 
hances the peak-to-background ratio, and does not create 
artifacts in the intensity distribution. 

Additionally, the Fermi surface maps in figures [^a, d), 
which span a portion of in-plane momentum space larger 
than in Fig. 2 of the main text, show that the periodicity 
of the electronic structure at the cleaved (001) and (101) 
surfaces, respectively, corresponds to the one expected at 
unreconstructed surfaces. 

Note that the Fermi surface map of Figure [^b), and 
the corresponding dispersions at hu = 47 eV {kz = 
3.45 A“^) shown in Fig. 3 of the main text, were not 
obtained at normal emission, but in the second Brillouin 
zone. Thus, the mere geometry of measurements imposes 
already extrinsic intensity asymmetries between the the 
left and right branches of the 2DEG electron bands. 
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FIG. 7. (Color online) (a) Raw Fermi surface map measured at hu = 47 eV on a cleaved, bulk insulating Ti02 anatase (001) surface. 
The red lines indicate the edges of the bulk Brillouin zones in the (001) plane of rio2- (b) Raw Fermi surface map in the A:^ooi) “ ^( 010)5 or 
(100) plane. The red lines are the bulk Brillouin zones containing rio2 and Fios. (c) Raw EDCs measured at hu = 47 eV (kz = 3.45 A“^) 
along the (010) direction on the (001) surface. They correspond to the energy-momentum map shown in Fig. 2(d) of the main text, (d) 
Raw Fermi surface map measured at hu = 47 eV on a cleaved, bulk insulating Ti02 anatase (101) surface. The red lines indicate the 
edges of the unreconstructed Brillouin zones in the (101) plane of roo2- (e) Raw Fermi surface map in the ky — kz plane measured on a 
cleaved, bulk insulating Ti02 anatase (101) surface. The red lines are the bulk Brillouin zones containing r2oo and Faoo- (f) Raw EDCs 
measured at hu = 47 eV (kz = 3.45 A“^) along the (010) direction on the (101) surface. They correspond to the energy-momentum map 
shown in Fig. 2(h) of the main text. 


Polarization dependence 

The polarization dependence of the t 2 g orbitals 
{dxy^dxz^dyz) forming the 2DEG at the (001) surface of 
SrTiOs was discussed in previous works [161SS] and can 
be directly applied to the (001) surface of anatase. To 
understand the photoemission selection rules, the geom¬ 
etry of the experimental setup is essential and shown in 
Fig. ga). We used light polarized horizontally (LH, in 
the x — z plane) and vertically (LV, along the ^-direction). 
The slits of the analyzer are aligned vertically along the 
^-direction and the sample was rotated around the x axis. 


In this geometry, the measurement plane containing the 
electrons ejected from the surface and entering the de¬ 
tector through the vertical slits, coincides with the yz 
plane. 

The measurements were conducted at a photon energy 
of hu = 47eV both at normal emission and tilted by 
26° around the x-axis. The corresponding positions in 
reciprocal space of bulk anatase, labelled A and B, are 
shown in figure [^b). 

To measure a non-zero intensity, the final state needs 
to be even with respect to the y — z measurement plane. 
As LV polarisation is even with respect to this plane. 















































< 010 > 


11 


incident 

light 



Anatase (001) 



^<010> (A ) 



) ..111111111111111111111111111111 p 

-0.4 -0.2 0.0 0.2 0.4 


k 


<oio> 


(A-^) 



(e) 


IiiiiIiiiiImiiIiimIiiiiIiiiiIiiiiIiimI 



I .I. I .I.I 


-0.4 -0.2 0.0 0.2 0.4 


k 


<oio> 


(A') 


FIG. 8 . (Color online) (a) Sketch of the experimental ARPES setup, (b) Schematic representation of the reciprocal space of bulk Ti02 
anatase over several Brillouin zones in the (100) plane, (c) Energy momentum map measured at hu = 47 eV with linear horizontal light 
polarisation (LH) at the surface of a tilted sample corresponding to position B in reciprocal space in (b). (d) Same as (c) but with linear 
vertical light polarisation (LV). (e) Same as (c) but at normal emission, position A in (b). 


only states with an orbital character even to the y — z 
plane, hence dyz states, can be detected using LV. The 
LH polarisation has even as well as odd components with 
respect to the y — z plane, and states of all t 2 g orbitals 
have non-zero matrix elements using this polarisation. 

The ARPES spectra in figures [^c, d) show the mea¬ 
sured intensities at the point B in reciprocal space for 
LH and LV polarizations, respectively. The almost zero 
intensity in LV polarisation demonstrates that the or¬ 
bital character is not dyz. To exclude one more orbital 
character we turn to the measured intensities at normal 
emission (point A) in Fig. [^e). To have a non-zero ma¬ 
trix element at the T point in normal emission, the final 


state has to be even with respect to both the x — z and 
y — z planes. Hence, as the LH polarization is even with 
respect to the x — z plane, only states with an orbital 
character even to the x — z plane {dxz) can be detected. 
The measured intensity at and close to normal emission 
in figure [^e), i.e. around /c = 0, is zero and thus, the 
bands are not of dxz character. 

Consequently, the orbital character of the bands is 
dxy, in agreement with the circular Fermi surfaces in fig¬ 
ure 2(b) of the main text and the non-degeneracy of the 
t 2 g orbitals in the bulk of anatase. We verified the ma¬ 
trix elements at a different photon energy to affirm that 
the matrix elements are due to orbital character and not 
photon energy. 





















